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A simple model for calculating the diffraction radiation characteristics 
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from an ultrarelativistic charged particle moving close to a tilted ideally con- 
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t/5 ', a superposition of the radiation fields for periodically spaced strips. The RDR 

characteristics have been calculated as a function of the number of grating 
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ducting strip is developed. Resonant diffraction radiation (RDR) is treated as 



elements, tilted angle, and initial particle energy. An analogy with both the 
resonant transition radiation in absorbing medium and the parametric X-ray 
radiation is noted. 
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acteristics of electromagnetic radiation from a relativistic particle moving parallel to and 
above diffraction grating (so-called Smith-Purcell effect (SPE)). The interest in this type of 



To date, a number of different approaches [1-4] have been employed to treat the char- 



radiation is related to both a possibility of using the SPE to generate intense radiation in 



the millimeter and submillimeter regions [5] and to implement it in non-destructive beam 
diagnostics [6] . In both cases we have to estimate the influence of such factors as transverse 
beam size, angular beam divergence, monochromaticity, etc. on the radiation characteristics. 
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In principle, the beam divergence effects in a plane that is parallel to the grating can 
be estimated using the results obtained in [4]. However, there is no simple algorithm to 
calculate the radiation characteristics from a particle passing above a tilted grating, i.e. for 
a non-zero angle between the particle trajectory and the main plane of the grating. 

One of the authors of the present paper has developed an approach based on description 
of SPE as a resonant diffraction radiation (RDR) [2] which is suitable for calculating the 
radiation characteristics for a grating consisting of a number of conducting strips spaced by 
vacuum gaps. Here we calculate the RDR characteristics of a particle whose trajectory is 
not parallel to the grating. Also, the influence of geometry on the RDR characteristics is 
studied. 

DIFFRACTION RADIATION (DR) FROM AN ULTRARELATIVISTIC 
PARTICLE FOR A TILTED STRIP 

To calculate the RDR characteristics for a tilted grating it is neccessary to know the 
diffraction radiation field from a charged particle moving close to a single tilted strip. 

The exact solution of Maxwell equations describing the radiation from a charged particle 
moving above an inclined semi- infinite ideally conducting screen has long been known [7]. 
Using the latter and the results reported in [8] , we can obtain an expression for the DR field 
strength for a tilted strip as a difference between the radiation fields of two semi-infinite 
planes, one restricted by edge 1, and the other — by edge 2 (see Fig.l) 

Hi strip tLiup tLi down \*-) 

For convinience, we express the impact parameter (the shortest distance between the 
particle trajectory and the plane edge) for edge 1 in the following way: 

a 
ai = h- -sine , (2) 

whereas for edge 2 



ai — h + — sin6 . (3) 

Here h is the spacing between the particle tragectory and the middle line of the strip, a 
is the strip width, and 0o is the strip tilt angle. Thus, taking into account the phase shift 
we have 

Estrv P =E D R (h - I sin o )e^- Edr (h + ^ sin e )e"^ (4) 

In (4) the DR field for a semi-infinite ideal screen is expressed through Edr {o-i}- The 
full phase shift (20) characterizes the phase difference between the waves being formed in 
the vicinity of edges 1 and 2 [9] ; and it can be derived from simple geometrical relations as 
a quantity which is propotional to the time difference between the wave propagation from 
edge 1 and edge 2 (see Fig.l) 
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20 = A *_A (5) 

where (3 is the particle velocity, and A is the DR wavelength. 

In this paper use is made of the system of units h — m — c — 1. 

For the extreme cases (Oq = 0° and 90°), the phase shift calculated according to (5) 
coincides with that obtained in [2]. 

Equation (4) is valid if the following conditions for the strip dimensions are fulfilled: 
7A >> b, 7 A << c (see Fig.l). Here 7 is the Lorentz - factor of the particle. 

We will consider the radiation characteristics in a coordinate system where the z axis 
is directed along the beam, the x axis is parallel to and the y axis is perpendicular to the 
strip edge. As shown in [8], for relativistic particles the radiation concentrates in the range 
of angles 

|e x | < 7" 1 , (6) 

if 7 >> 1. 

In this approximation the following expression for the DR yield with wavelength A will 
be true: 



Edr (x + Ax) =Edr (x) exp 



2TrXr 'i + f»e* 



(7) 
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where Edr (x) ~ exp j — ^7- y 1 + 7 2 © 2 }, anc ^ Uc = 2o~ * s ^^ e ^^ characteristic energy. 
For further calculations we shall use a more symmetrical expression instead of (4) 

Estri P =EDR (h) [exp(a + i(f>) - exp(-a - i(f>)] , (8) 
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a = ^ ^ ° j y/l + 7 2 © 2 (9) 

From (8) we derive the following fomula for the DR spectral-angular density for the strip: 

d 2 W strip d 2 W DR 



duodQ, duodVt 



■F str , (10) 



d 2 W, 



D " 4n 2 \ Edr \ 2 ,F str = 4 (sinh 2 a + sin 2 (f)) (11) 



duidQ 

The expressions obtained are quite similar to the formulas for the spectral - angular 
distribution of transition radiation (TR) from a foil (see, e.g., [10]). In the case in question, 
F str characterizes the DR field interference from the two strip edges; whereas, for TR, the 
same multiplier characterizes the TR field interference from the input and output surfaces 
of the foil. 

Earlier [9] , the spectral-angular density of DR for a semi-infinite screen has been obtained 
using an ultrarelativistic approximation. It was shown that the DR is concentrated in the 
vicinity of the plane that is perpendicular to the screen (Q x ~ 7 _1 ) 



S a = 4^ eX P(-^V /l+ ^ X 
X {02[l +C os(0 J/ -0 o )](l-COS0o) + 

{T 2 + ©x)[l - cos(e w - O )](1 + cos ©o)} x 

x {( 7 - 2 + e»){[oos(e s -e )-^] 2 + 

+ ( 7 - 2 + © 2 )sin 2 e }}- 1 (12) 

In (12), we omitted the terms smaller than ~ 7 -2 in the numerator and denominator the 
terms. 



Now we shall consider the forward diffraction radiation (FDR), i.e., for the angles Q y ~ 
7 _1 << f . Using this approximation instead of (12) we have 

d 2 W DR a ( u r ; 2r ^ 2 \ 7~ 2 + 20 2 

ex P —V 1 + 7 2 © 2 7—0 , 75^7—0 I Q2 I ^ ( 13 ) 



dudtt 47T 2 ^ W C V ' X J ( 7 -2 + e2)( 7 -2 + 62+02) 

As was noted in [8], in the angular distribution of DR at Q y = Q x = there is a maximum 
whose value is proportional to 7 2 



d 2 W DR (Q x = 0, e„ = 0) _ a 
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c^n A^ exp {-uJ (14) 



For G y >> 7 *, using (12) one can obtain a simpler expression: 
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^m = -S-exp (-^,/l+7 2 62 



x { 7 - 2 (l + cos O )[1- cos(Q y -Q )] + 

+ 20 2 [1 - cos O cos(0 y - 6 )]} x 
x {(T 2 + ©') sin 2 ^ sin 2 (e - ^) T 1 (15) 

As follows from (12), for the mirror reflection angle (Q y = 20o) the angular distribution 
has another maximum which is much the same as Eq. (14) and can be identified as "backward 
diffraction radiation" (BDR) in analogy to the process of transition radiation. The fact that 
the FDR and BDR intensities coincide in all the frequency range results from using an ideally 
conducting screen approximation. It is obvious that, however, the latter approximation is 
not valid for uo c > uj p , where uo p is the plasmon energy of the screen material. 

In the range of angles |20 o — Q y \ ~ 7 -1 , i.e., in the vicinity of the mirror reflection 
direction, the BDR spectral-angular density also has the form (13); however, in this case 
the angle Q y is measured from the mirror reflection direction (see [8]). 

Thus, expression (15) is valid when the following conditions are fulfilled: 

0, »7~ 1 , |20 o -0,|»7 _1 (16) 

In this case, as follows from (15), the DR density is about 7 2 times smaller. 



Figure 2 shows the dependence of the DR yield on the strip tilt angle O (so-called 
orientation dependence) for the fixed observation angle Q y = 4.5° in the reflection plane 
(Q x — 0)- The calculations have been carried out using the formulas (10) — (12) for 7 = 1000, 
A = 0.4/xm, and 01 = 0.1mm. One can notice that for 7A > a sin Oo the DR yield is strongly 
suppressed and the characteristic angular width of the dependence is significantly higher 
than ~ 7 _1 which is typical for DR from a semi-infinite screen (see Fig. 2d). The DR 
intensity at the orientation dependence maximum is also suppressed when 7A > asin@o- 

RESONANT DIFFRACTION RADIATION FROM A TILTED GRATING 

Let us consider a grating consisting of N strips, of width o and period d tilted at angle 
0o to the electron momentum (see Fig. 3). The impact parameter (the distance between the 
first strip center and the electron trajectory (1)) we denote as h. 

The radiation field being formed near strip 1 of the grating coincides with expression (8) 

Ei=Estr (h) (17) 

whereas the next strip field differs from (17) by both the phase 0o an d the decay factor 
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Ei=E\ exp(-a - i(f) ) =E s tr (h) exp(-a - i(f) ), (18) 



which can be determined in analogy to (5) and (9) 

2nd [cos^ - O ) - ^^ 



A 
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(19) 



One can write the k th strip field in the same way 

Ek=Estrexp[-(k-l)(a + i<f)o)] (21) 

The resulting field of the iV strip grating is expressed through the sum of iV terms 
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Egr=Ei + i 2 + • • • + E k =Estr ]T exp[-(fc - l)(a + i</> )] (22) 

fc=i 

Having calculated the squared modulus of expression (22) we obtain the following ex- 
pression for the DR spectral-angular density for the entire grating: 

d 2 W GR d 2 W str d 2 W DR 



dujdVt dudVt ' dudVt 



where 



1 N 



N 



J2 exp[-(fc - l)(a + i(j>o)] 
fc=i 



l-C N 2 



C 



(24) 



Here C = exp(— a — z0o) 

After simple mathematical transformations, expression (24) may be written in the fol- 
lowing manner: 

One should notice that the structure of formula (25) is identical to that of a similar ex- 
pression for the resonant transition radiation from N layers taking into account the radiation 
absorption in every layer. 

First we shall consider a particular case corresponding to the Smith-Purcell geometry 
(Go = 0). It is obvious that the decay factor «o determied by expression (20) is equal to 
zero, therefore, expression (25) can be rewritten in a well-known form 

sin 2 (^) , x 

F " = —hU ( 26 ) 

sin (*f) 

When N — *■ oo, expression (26) transforms into an ordinary 5- function 

F N = 2nN8((f)o - 2kn) (27) 

where k is the diffraction order. 

The presence of 5- function is an indication of the existance of the monochromatic maxima 
in the RDR spectrum. However, the use of the S— function for real gratings where the 
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number of elements is limited is not always justified. Therefore, we shall further use the 
exact formulas (25) and (26). 

Figure 4 depicts the RDR spectral distribution for the Smith-Purcell geometry. The 
calculation has been carried out for a = | when the intensity reaches its maximum value 
[2]. The peak position in the spectrum is determined by the phase relation (the resonance 
condition) 

0o = 2kir, (28) 

which leads to a well-known formula of Smith-Purcell 

d (cos©,, — 4 ) 

Afc = -i J—H (29) 

k 

As is seen from the figure, the even orders are absent. This is explained by the influence 
of the F str factor that is equal to zero at 

(j) = mil, (30) 

where m is the integer (see formula (11)). 

Substituting (29) in (5), we can write expresion (28) in the following form: 

d 

k = —m = 2m (31) 

a 

Thus, in the case under consideration (a — |) the even diffraction orders are forbidden. 
For illustration, figure 5 shows the dependence of F str on the photon energy which was 
calculated for the same conditions as in Fig. 4. 

Figure 6 shows the maximum DR yield dependence on the ratio | for four radiation 
orders . 

As was shown above, for the first diffraction order the intensity reaches its maximum 
value at § = 0.5, whereas for the second order — at § = 0.25 and 0.75. 

Let us consider the case of a tilted grating. One should note that in the DR spectrum 
the quasimonochromatic peaks can be observed at small tilt angles of the grating. For the 



tilted single strip the spectrum calculated according to (10) is shown in Fig. 7. The strip 
parameters and geometry are indicated in the figure caption. As in the case of a semi-infinite 
screen one can observe an exponentially decreasing spectrum. 

Figure 8 depicts the RDR spectrum from a tilted grating calculated at Go = 1.9°. 
Unlike the spectrum from a single strip calculated for the same initial conditions, here one 
can observe the first order quasimonochromatic radiation maximum with a finite full width 
at half maximum (FWHM) together with the continuous background. 

The position of the quasimonochromatic maxima in the RDR spectrum is determined 
by the resonance condition (28) where the <pQ phase is taken according to formula (19). We 
will illustrate this fact in the following way. Let us rewrite expression (25) in the form 

l-2exp{-(iV-l)ao}cos{(iV- 1)<M + exp{-2(7V - l)a } 



1 — 2 exp{— a } cos 0o + exp{— 2a } 
In extreme case, when N — ► oo, instead of (32) we have 

1 



(32) 



(33) 



1 - 2e~ a ° cos 0o + e- 2a o 

It is apparent that the expression obtained reaches its maximum value when the following 
conditions are fulfilled: 

<p = 2kn (cos0 o = l) (34) 

In this case: 

F »=(id^ (35) 

At small values, ao << 1, that correspond to small tilt angles of the grating, 6o << 1, 
from (35) we have 

Foo = \ (36) 

"o 

Thus, the pronounced quasimonochromatic maxima in the DR spectrum can be observed 

at the sliding incidence angles of the particle beam with respect to the grating. For the 



angle Oo 7^ 0, instead of Smith-Purcell conditions we have the following relation between 
the quasimonochromatic maximum position, period d, grating tilt angle 0o and observation 
angle Q y : 

d 



Xk 



008(9, - e„) - ss^tt 



(37) 



k 

where k is the integer. 

Figure 9a shows the dependence of different - order maximum positions on the target 
orientation angle at the observation angle Q y = 4.5°. 

One can notice that for negative values of Oo (i.e. for the geometry where the beam 
"reflected" by the grating is directed to the opposite side from the detector) the spectral 
maxima are shifted in the soft part with respect to the SPE spectrum (Bo = 0). When the 
tilt angle increases, the spectral maxima are shifted to the high energy part. 

A similar peak shift of the parametric X-ray radiation has been registered in the exper- 
iment [11] with rotating a crystallic target. Actually, for both the RDR and PXR the peak 
position in the spectrum is determined by the resonance condition only and does not depend 
on the radiation mechanism. 

The experiments on the PXR studies often measure the so-called theta-scan, i.e., the 
dependence of the radiation yield for a fixed observation angle and electron energy on the 
target orientation angle. Figure 9b presents this dependence calculated for the RDR. In this 
case the grating tilt angle Go is varied with respect to the electron beam. 

A recent experiment [12] measured the same dependence for the Smith-Purcell effect for 
the grating made as a periodically deformed continuous surface. The authors of the quoted 
paper obtained the dependence with clear maxima. This dependence shape is quite close to 
the one presented in Fig. 9b. 

Figure 10a presents the dependence of the full width at half maximum (AE) on the 
number of grating elements for the first diffraction order and tilt angle 0o = 1° (solid line). 

It also shows a similar dependence for the SPE which is well approximated by the -^ 
dependence (dotted line). As follows from the figure, for G = 1° the resulting curve is well 
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approximated by the ^ formula where C\ = 1.8. 

When the number of periods increases, the RDR intensity also increases reaching the 
0.95/oc level for N = 70 (see Fig.lOb). 

Let us estimate the effective grating length (the number of periods) with which the 
passing particle field interacts 

N *ft = T^TV ( 38 ) 

d sin O 

In the case considered N e ff = 60, which is quite close to the grating length providing a 
95% intensity level. 

When the particle moves close to the grating of a limited length (with the number of 
periods N), one can derive a similar caracteristics for the Lorentz- factor 

d sin ©o , 

left = N ^ — ( 39 ) 

Figure 11 shows the dependence of I max on the particle energy. As follows from the figure 
the simple estimation (39) is a good characteristics of the RDR process too. 

SUMMARY 

A simple model for calculating the RDR characteristics from a tilted grating has been 
suggested. It has been shown that quasimonochromatic maxima appear in the RDR spec- 
trum, and their characteristics (full width, intensity) are determined primarily by the angle 
between the grating plane and electron pulse. 

Let us estimate the initial beam divergence effect for the SP geometry. It was noticed 
in [4] that for the angles 0# << 1 (0# is the angle between the electron momentum and 
grating axis in the horizontal plane) the peak shift can be described in the following way: 

cos0# k \ 2 J 

Thus, 
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therefore, the peak broadening related to the beam devergence in horizontal plane can 

be neglected if 

2 1 

^ « -, (42) 

2 AT' v ; 

which is generally fulfilled. 

In order to obtain a formula, analogous to (41) and characterizing the beam broadening 
due to the divergence in the vertical plane, we shall take into account that the observation 
angle 

@jz — @o — @d = const, (43) 

where the angle 0£> is measured from the grating plane. 

Let us denote the angle 0o as 0y (0y is the angle between the electron pulse and grating 
axis in the vertical plane). Then for 0y << 1 we have 



A ' 02 

\ - d 
Xk ~k 



D — ■ 



COS - -' 



(44) 



It follows from (44) that the SPE peak broadening is determined by the observation 
angle 0£> 



AX V 0? 

A " 4sin 2 % 



(45) 

2 

For the small observation angles 0£> the peak broadening ^^ can be significant and 
exceed the "natural" peak width jj. Thus, changing the SPR peak width, for the small 
angles 0d, in principle, one can determine the vertical beam divergence. 

From (45) one can estimate the sensitivity of the method: 

2 sin % 

e, > -J- (46) 



Since the latter expression does not depend on the wavelength, investigation of the 
shape line can be carried out in the optical region. For example, for 0£> = 2° and N = 100 
0y ~ 3.5 ■ 10~ 3 . The estimation obtained does not depend on 7 (if 7 >> 1). 
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The suggested technique for the beam divergence determination can be used in accelera- 
tors with 7 < 100, because, in this case, the well-known methods based on either transition 
radiation or synchrotron radiation exhibit a sensitivity controlled by the characteristic ra- 
diation angle 7~ 1 . 

The dependence of the peak position in the RDR spectrum on the grating tilt angle 0o 
can be used to determine the electron bunch length l e . In [6,13] the authors suggested to 
measure the coherent SPR yield at different observation angles. In the wavelength region of 
A ~ l e one will observe conversion from an ordinary SPR to coherent SPR, i.e., the radiation 
intensity will be changed by approximately N e times (iV e is the number of electrons in a 
bunch) . 

Figure 12a shows the dependence of the RDR wavelength on the grating tilt angle for 
the fixed observation angle Q y , and a similar dependence for the peak width is presented in 
Fig. 12b. The calculations have been carried out for the following conditions: 

7 = 100, Q y = 4.5°, d = 20mm, a = 10mm. 

As follows from the figure, in the wavelength region A = 0.15 -f- 0.8mm(O = 8° -j- 35°) 
a practically constant intensity of the ordinary RDR is observed. If we measure the RDR 
intensity for the electron bunch with the length l e ~ 0.5mm, then for the grating tilt 
angles indicated one could research the conversion to the "coherent RDR mode" in detail, 
which would allow one to determine both the everage bunch length and the profile electron 
distribution in a bunch. The suggested technique is related to the grating rotation for a 
fixed detector position, whereas in [6,13] it was sugested to move the detector, which is not 
always convinient (in particular, if the detecting system contains a monochromator). 

Depending on the bunch length l e , the region of the investegated wavelengths can be 
easily changed by choosing a proper grating period and observation angle. 
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FIGURE CAPTIONS 

Fig.l. Geometry of diffraction radiation (DR) from a single strip; h is the impact pa- 
rameter, Q y is the observation angle, Gq is the strip tilt angle, and a is the strip width. 
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Fig. 2. Orientation dependence for a single strip with the widths a=0.1mm (a), 1mm (b), 
10mm (c), and a semi-infinite plane (d). 

Fig. 3. Geometry of diffraction radiation from a tilted grating; h is the impact parameter, 
a is the sprip width, d is the grating period, and G>o is the grating tilt angle. 

Fig. 4. Smith-Purcell effect spectrum. The initial conditions used are Q y = 4.5°, Q x = 0, 
0o = 0, a = 0.2mm, d = 0.4mm, h = 0.1mm, 7 = 1000, N = 50; k is the diffraction order. 

Fig. 5. Dependence of F STR on energy. 

Fig. 6. Dependence of the first (1), second (2), third (3) and fourth (4) order maximum 
intensity on the ratio |. The initial conditions used are Q y = 4.5°, O^ = 0, G>o = 0, 
h = 0.1mm, 7 = 1000, N = 50, A = 0.4/zm. 

Fig. 7. Diffraction radiation spectrum from a tilted strip. The initial conditions used are 
Q y = 4.5°, 9^ = 0, O = 1.9°, a = 0.2mm, 01 = 0.1mm (see Fig.l.), 7 = 1000. 

Fig. 8. Diffraction radiation spectrum from a tilted grating. The initial conditions used 
are & y = 4.5°, 9,, = 0, 6 = 1.9°, a = 0.2mm, h = 0.1mm (see Fig.3.), 7 = 1000. 

Fig. 9. a: Dependence of the first (1), second (2), third (3) and fourth (4) order peak 
positions (Q y = 4.5°, d = 0.4mm) on the tilt angle. 

b: Dependence of the first DR order yield (A = 1233nm) on the tilt angle. 

Fig.10. a: Dependence of the first order FWHM for SPE (0 O = 0) (solid line), O = 1° 
(dash-dotted line) on the number of elements and dependence of the -^ (dashed line). The 
initial conditions used are Q y = 4.5°, Q x = 0,a = 0.2mm, h = 0.1mm, d = 0.4mm (see 
Fig.l.), 7 = 1000. 

b: Dependence the first order maximum intensity on the number of elements 
calculated for the same initial conditions. 
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Fig.ll. a: Dependence of the first order FWHM (dash-dotted line) on the N e ff and 
dependence of j^— (solid line). The initial conditions used are 6^ = 4.5°, Q x = 0, ©o = 1°, 
a = 0.2mm, h = 0.1mm, d = 0.4mm (see Fig.l.), 7 = 1000. 

b: Dependence of the first order maximum intensity on the N e ff calculated for 
the same initial conditions. 

Fig. 12. a: Dependences of the RDR peak intensity on the grating tilt angle (radiation 
wavelength) for 7 = 100, Q y = 4.5°, d = 20mm, a = 10mm. 

b: Dependence of the FWHM on the grating tilt angle (radiation wavelength) 
for the same initial conditions. 
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